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Abstract: In this note we investigate what is the best LP-norm in order to describe the relation 
between the evolution of the state of a bilinear quantum system with the L^'-norm of the external 
field. Although has a structure more easy to handle, the norm is more suitable for this 
purpose. Indeed for every p > 1 it is possible to steer, with arbitrary precision, a generic bilinear 
quantum system from any eigenstate of the free Hamiltonian to any other with a control of 
arbitrary small norm. Explicit optimal costs for the norm are computed on an example. 
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1. INTRODUCTION 

1.1 Physical context 

The state of a quantum system evolving on a Riemannian 
manifold f], with associated measure fi, is described by 
its wave function, that is, a point in the unit sphere of 
L^(rj, C). A system with wave function is in a subset w 



zero. When the control u is nonzero, and provided suitable 
regularity hypotheses, the energy evolves as 



^ = 2u{t)^ (^J^ [(A + V)iJ\ Wip dfi 



(2) 



of O with probability \^\ d/i. 

When submitted to an excitation by an external field (e.g. 
a laser) the time evolution of the wave function is governed 
by the bilinear Schrodinger equation 

i^ = -AV- + V{x)^P{x, t) + u{t)W{x)^{x, t), (1) 
at 

where V,W : VL ^ Tl are real functions describing 
respectively the physical properties of the uncontrolled 
system and the external field, and u : R — >■ R is a real 
function of the time representing the intensity of the latter. 



1.2 Energy for a quantum system 

Physically, the energy of a quantum system (1) with wave 
function ^ is £'(^) = [(~^ + V' d/i. The energy 
is therefore constant in time whenever the control u is 



Note that the time derivative of the energy E at time t 
depends on the value u{t) of the intensity of the external 
field and on the wave function ip{t). Therefore a natural 
question is to find an a priori relation between the time 
evolution of the energy of system (1) and properties 
of the external field represented by u. In particular we 
address the problem of finding a bound on the energy after 
the action of an external field. Namely, given an initial 
condition ij}^ and a control m : [0, T] — >■ R denoting with ip 
the solution of (1) with = V'o we look for bounds on 
the energy E{T) in terms of the norm of u 



l|w|Uf(0,T) 




u{t)\Pdt 
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(for some suitable p > 0) without computing explicitly the 
solution t e [0,r] i-> il){t). 

Many previous works adressed the problem of the optimal 
control of the system (1) for costs involving the norm 
of the control (see for instance Dahleh et al. (1990) or 
Grivopoulos and Bamieh (2008)). The main reason for the 
choice of the norm is the fact that the natural Hilbert 
structure of spaces allows to use the powerful tools of 
Hilbert optimization. It is sometimes believed that there 
is a natural relation of the LP' norm of u and the energy 
of the systems. This note presents a priori bounds on the 
L^-norm of the control and shows that, in general, the 
L^-norm provides more informations on the evolution of 
the systemthan other L^'-norms for p > \. Most of the 



material presented below is classical for finite dimensional 
conservative bilinear systems. The contribution of the 
present note is to treat in a rigorous and unified way the 
case of both finite and infinite dimensional systems. 

1.3 Content of the paper 

The first part of the paper (Section 2) presents some 
theoretical tools for bilinear quantum systems. The ques- 
tion of estimation of the energy is reformulated in terms 
of a problem of optimal control. Some relations between 
the variation of energy of a quantum system and the 
norm of the external field are given in Section 3. Finally, 
some explicit computations are presented on an example 
in Section 4. 

2. INFINITE DIMENSIONAL QUANTUM SYSTEMS 
2.1 Abstract framework 

We reformulate the problem (1) in a more abstract frame- 
work. This will allow us to treat examples slightly more 
general than (1), for instance, the example in (Boussai'd 
et al., 2011b, Section III. A). In a separable Hilbert space 
H endowed with norm || • || and Hilbert product (•,•), we 
consider the evolution problem 

^ = (A + H(t)B)^ (3) 

where {A, B) satisfies the following assumption. 
Assumption 1. {A,B) is a pair of linear operators such 
that 

(1) A is skew-adjoint and has purely discrete spectrum 
(— iAfc)fcgNi the sequence (Afc)/.gN is positive non- 
decreasing and accumulates at -|-c»; 

(2) B : H ^ H \s skew-adjoint and bounded. 

In the rest of our study, we denote by {(j)k)k£'H an Hilbert 
basis of H such that = — iAfct^^ for every k in N. 

We denote by D{A + uB) the domain where A + uB is 
skew-adjoint. 

From Assumption 1.2 together with Kato-Rellich Theo- 
rem, we deduce that A -\- uB is skew-adjoint with domain 
D{A). Moreover, for every constant u in R, \{A + uB) is 
bounded from below. 

Hence, for every initial condition ipQ in H, for every u 
piecewise constant, u : t ^ J2j '^jXitj,tj+i){t), with — 
to ti < . . . < tAT+i and uo, . . . ,U]\r in R, one can define 
the solution t Tfipo of (3) by 

o g(tj-i-tj-2)(A+«,_2-B) o . . . o e*"(^+"°^)Vo, 

for t € [tj-i,tj). For a control u in i^(R) we define the 
solution using the following continuity result. 

Proposition 1. Let u and (w„)„gN be in L^(R). If for 
every tin R u„(r)dT converges to M(r)dT as n tends 
to infinity, then, for every t in R and every -tpQ in H, 
{Tf"ipo)ne'!Si converges to T^ipo as n tends to infinity. 



controllability of (1) in is hopeless in general. With the 
exception of some very particular examples where Q is one 
dimensional, (Beauchard and Laurent (2010)), no descrip- 
tion of the attainable set is known. Therefore, we often 
consider the weaker notion of approximate controllability: 

Definition 2. Let {A, B) satisfy Assumption 1. The system 
{A^B) is approximately controllable if, for every '0Oj'0i 
in the unit Hilbert sphere, for every e > 0, there exists 
Ue : [0,Te] ^ R such that ||T^J?/>o - V^iH < e. 

Various methods have been used to give criterion of ap- 
proximate controllability of system (1). Nersesyan (2010), 
Beauchard and Mirrahimi (2009) and Mirrahimi (2009) 
rely on a Lyapunov approach. Boscain et al. (2012) adopt 
a more geometrical point of view, centered on the notion 
of non-degenerate (or non-resonant) transitions. 

Definitions. Let {A,B) satisfy Assumption 1. A pair 
(j, k) of integers is a non- degenerate transition of (A, B) 
if (i) {(j)j,B(j)k) 7^ and (ii) for every {l,m) in N^, \Xj — 
Afel = |Ai - Ami implies {j, k) = {I, m) or {(jyi, Bcpyn) = or 
{j, fc}n{Z,m} = 0. 

Definition 4- Let [A^B) satisfy Assumption 1. A subset 
5* of is a non- degenerate chain of connectedness of 
{A, B) if (i) for every (j, k) in S, (j, k) is a non-degenerate 
transition of {A,B) and (ii) for every ra,rb in N, there 
exists a finite sequence = rp, ri, . . . , rp = rt, in N such 
that, for every j < p — 1, {rj, rj+i) belongs to S. 

The following sufficient criterion for approximate control- 
lability is the central result of Boscain et al. (2012). 

Proposition 5. Let {A^B) satisfy Assumption 1. If {A,B) 
admits a non-degenerate chain of connectedness, then 
(A, B) is approximately controllable. 

As proved by Mason and Sigalotti (2010) (see also Privat 
and Sigalotti (2010)), a system {A,B) satisfying Assump- 
tion 1 generically admits a non-degenerate chain of con- 
nectedness. Hence, approximate controllability is a generic 
property for systems of the type (3) (see also Nersesyan 
(2010)). 

In Boscain et al. (2012) an extensive (and in some case 
implicit) use of averaging results has been made. The 
following result is a generalization of the Rotating Wave 
Approximation to infinite dimensional systems and can be 
found in Chambrion (2011). 

Proposition 6. Let [A, B) satisfy Assumption 1 and (j, k) 
be a non-degenerate transition of {A,B). Define T = 
27r/|A, - Afel and = {{l,m) e N2|(0,,B0,„) ^ 
and |A,-A,„| e (N \ {1})|A, - Afc| and {j,fc}n{/,m} ^ 
0}. Consider a T-periodic function u* : R — > R satisfying 

/o^w*(t)e'(^^-^'=)*dt ^ and /J' (t)e'(^'-^")*dt = for 
every {I, m) in Af and let 



ttT 



2|6i,2| /o ^t*(T)ei(^l-^2)rd^- 



^.^ Controllability results 

Considerable efforts have been made to study the control- 
lability of (1). It is known (see Turinici (2000)) that exact 



Then there exist a sequence (r^)„gN such that T* € 
{nT* - T,nT* + T) and |(0a;, T^7"0j)| tends to one as 
71 tends to infinity. 



2.3 Formulation of an optimal control problem 

Let ( A, B) satisfy Assumption 1 and admit a non- 
degenerate chain of connectedness. For every r > 0, for 
every j, fc in N and e > we define A^{j, k) as the set of 
functions u : [0, T„] ^ R in L^{[0,Tu]) L''([0,T„]) such 
- < e. We consider the quantity 



that ||T^^(/)j 



sup I inf ||m||l'-(o.t„) 

e>0 \neA'^{],k) 



This quantity is the infimum of the L''-norm of a control 
achieving approximate controllabihty. It clearly satisfies 
the triangle inequality. Next proposition states that Cr is 
a distance on the space of eigenlevels only when r = 1. Its 
proof is given in Section 3. 

Proposition 7. Ci is a distance on the set {(j>j,j G N}. For 
r > 1, Cr is equal to zero on the set {(f>j,j G N}. 

2.4 Weakly- coupled systems 

Definition 1. Let fc be a positive number and let {A,B) 
satisfy Assumption 1. Then {A,B) is k weakly- coupled if 
for every m € R, D{\A + uiB\'^/'^) = D(\A\''/'^) and there 
exists a constant such that, for every ij} in Z?(|A|'^), 

|3?(|A|V,BV')I <C(A.B)|(|A|V,V')|- 

The notion of weakly-coupled systems is closely related to 
the growth of the |A|'^/^-norm For fc = 1, this 

quantity is the expected value of the energy of the system. 

Proposition 8. (Boussai'd et al., 2011b, Proposition 2) 
Let {A,B) be fc-weakly-coupled. Then, for every t/^o S 
D{\A\''/^), K > 0,T >0, and u in L^{[0,oo)) for which 
\\u\\l. < K, one has |lT^(^o)|ifc/2 < e=(^^^)^|lV'o||fe/2. 

For every N in N, we define Cm the linear space spanned 
by (/>!, 02, . . . , 0Af and tt^t : H ^ H, the orthogonal 
projection onto Cn- 

N 
k=l 

The compressions of order N oi A and B are the finite rank 
operators A^^^ = ttnA^Cn S^^^ = t^nB^Cn- The 
Galerkin approximation of (3) at order N is the infinite 
dimensional system 

— x = + (4) 

at 

Since Cn is invariant by (4), one may also consider (4) as 
a finite-dimensional system, whose propagator is denoted 
by 

Proposition 9. (Boussai'd et al., 2011b, Proposition 4) Let 
fc and s be non-negative numbers with < s < fc. 
Let {A, B) be fc weakly-coupled Assume that there exists 
d > 0, < r < fc such that \\BiIj\\ < d\\ip\\r/2 for every ip 
in D(|A|''/2). Then for every e > 0, A' > 0, n e N, and 
itpj)i<j<n in £'(1^1''/^)" there exists N € N such that for 
every piecewise constant function u 

\\u\\li < If ^ ||Tr(V,) -^(V)(i,0)7rAr^,!U/2 < e, 
for every t > and j — I, . . . ,n. 

Remark 1. An interesting feature of Propositions 8 and 
9 is the fact that the bound of the norm of 



the solution of (3) or the bound on the error between 
the infinite dimensional system and its finite dimensional 
approximation only depend on the norm of the control, 
not on the time. 



3. PROOF OF PROPOSITION 7 

3.1 Lower hounds for the norm 

Proposition 10. Let (A, B) satisfy Assumption 1. For ev- 
ery j, fc in N such that B(j)j 7^ 0, for every locally integrable 
u : [0,T] R, 



ll"IUi(0,T) > 



Proof: Let j, fc in N. For every locally integrable u : 
[0,T] R, define y : u <^ e-^'TJ'^fc. For al- 

most every t, y is differentiable with respect to t and 
dy 

dt 

^ < \u{t)\\\i^ipj\ 



□ 



^ = u{t)e ^^Be'^^y. In particular. 



dt 

This concludes the proof of Proposition 10. 

A consequence of Proposition 10 is that Ci(0j,0fc) is 
bounded away from zero as soon as j 7^ fc. 

3.2 Upper bound for tfie norm 

In order to give an upper bound for Ci(0j, (j)k) when (j, fc) 
is a non-degenerate transition of {A, B), we come back to 
Proposition 6. In the case where J\f is finite, Boscain et al. 
(2012) give an explicit construction of a piecewise constant 
u* with value in [0, 1] and satisfying the assumptions of 
Proposition 6 such that, for every n in N, 

< 



Li(o,r*) 

More details about the choice of u* when {A, B) is weakly- 
coupled are given by (Boussai'd et al., 2011a, Section III- 
C). Let us just mention that the choice w* : t H> cos(|Aj — 
A|t) guarantees 



< 



1(0. 



for every n in N. 

In any case, {A, B) being weakly-coupled or not, this 
guarantees that, for every j, fc in N, Ci{(f>j,(t)k) < +00. 

The fact that Ci is symmetric is a consequence of the so- 
called time reversibility of (3): the propagator associated 
with {—A, —B) is the adjoint of the backward propagator 
associated with {A,B). This proves that Ci is a distance 
on the set {(f>j,j E N}. 

3.3 U" norms with r > 1 

Proposition 11. Let {A^B) satisfy Assumption 1 and ad- 
mit a non-degenerate chain of connectedness. If r > 1, 
then, for every j, fc in N, Cr{4>j, (pk) =0. 



Proof: It is enough to consider (j, k) in a non-degenerate 
chain of connectedness of {A,B). The resuh is then a 
consequence of Proposition 6, since, for every n in N, 



1 



L'-(0,T*) 



< 



n' 



nT*+T 

\u*{t)Ydt 







1 /•("[t^I+i)?' 

< — / \u*{t)Y>it 

k) 

< — ( ^ + 2) l\u*{t)Ydt, 



n' 

n fT* 
n' \ Y 



which tends to zero as n tends to infinity. 

4. ROTATION OF A PLANAR MOLECULE 



□ 



In this Section, wc apply our rcsuhs to the wcU studied 
example of the rotation of a planar molecule (see, for 
instance, Salomon and Turinici (2005); Boscain et al. 
(2009, 2012)). 

J^.l Presentation of the model 

We consider a linear molecule with fixed length and center 
of mass. We assume that the molecule is constrained to 
stay in a fixed plane and that its only degree of freedom 
is the rotation, in the plane, around its center of mass. 
The state of the system at time t is described by a point 
9 ^ ip{t,e) of L'^{n,C) where il = R/27rZ is the one 
dimensional torus. The Schrodinger equation writes 

i^(t, 0) = -Aiit, 9) + u{t) cos(0)V(i, 0), (5) 

where A is the Laplace-Beltrami operator on Vl. The 
self-adjoint operator —A has purely discrete spectrum 
{/c^,/c G N}. All its eigenvalues are double but zero 
which is simple. The eigenvalue zero is associated with 
the constant functions. The eigenvalue k"^ for /c > is 
associated with the two eigenfunctions 9 i-> cos{k9) 

and 9 -^sm{k9). The Hilbert space H = L'^{VL,C) 

splits in two subspaces and Ho, the spaces of even and 
odd functions of H respectively. The spaces and Ho 
arc invariant under the dynamics of (5), hence no global 
controllability is to be expected in H . 

We focus on the space Ho- The restriction A of iA to 
Ho is skew adjoint, with simple eigenvalues (— ifc^)fcgN 
associated with the eigenvectors 



s\n{k9) 



feeN 



The restriction _B of -0 i-> — icos(^^)■^/; to Ho is skew- 
adjoint and bounded. The pair {A, B) satisfies Assumption 
1 and is weakly-coupled (see (Boussai'd et al., 2011b, 
Section III.C)). 

The Galcrkin approximations of A and B of order N are 
/i ••• \ 



4i 





iV^i / 



and 



/ 1/2 
1/2 1/2 
'■• 



\ 







V 



1/2 1/2 
1/2 0/ 



4-2 Computation of Ci 



2) 



Our aim is to compute the minimal norm needed to 
approximately transfer the wave function from the first 
eigenspace to the second one. Precisely, we will prove 

Proposition 12. Ci(0i,02) = tt. 

Proof: The transition (1,2) is non-degenerate. Proposi- 
tion 6 applies with A/" = 0. For every 77 in (0, 1), we define 
, the 27r/3 periodic function defined by 

w''(x) = 1 for < X < 77 
w''(x) =0 for r/ < x < 27r/3. 

Proposition 6 states that, defining 

T* = 

|sin(f)|' 

there exists a sequence (r^)„gN such that T* € {nT* — 



27r/3,nr* + 27r/3) and |((?!)2,T^_ 
tends to infinity. 

One computes, for every n in N, 



< 



tends to one as 



37r 



sm 



ii(0,T*) 

This last quantity tends to tt as 77 tends to zero. This proves 
that Ci{(f>i, (j)2) < TT. 

For every locally integrablc control u, wc define yi : 
t {(j)i,Tf(j)i) and y2 ■■ t ^ {(j)2,T^(j)i). The function 
i I— >■ |yi(i)P is absolutely continuous, and, for almost every 
t, 

^^\y^it)\' = 2^iy[it)y^{t)) 



Hence, 



or 



= um{y2{t)m{t)). 
<\y2myi{t)Mt)\, 

< Ht)\ 



Integrating between and T, one gets 



arctan 



and, provided ||u||li(o,t) < ""j 
I 1 

— 1 < tan 

"ll-Li(0,T) 



V \yiiTW 

Finally, \yi{T)\ > cos 



\y2{T)\<Vi-\yi{tW< 




If (Wn)neN is & Sequence of locally integrable functions 
and (T„)„gN is a sequence of positive numbers such that 
\{(f)2,Tlp^ tends to one, then liniinf„ ||u„||^1(-o^t„) ^ 
TT, hence Ci{(j)i, (l>2) > ti"- This concludes the proof of 
Proposition 12. □ 

5. CONCLUSION AND PERSPECTIVES 

5.1 Conclusion 

We introduced an optimal control problem associated with 
a bilinear quantum system. For p > I, we have given some 
estimates of the norm of the control needed to steer 
the system from an eigenstate of the free Hamiltonian 
to another. In particular, for generic bilinear quantum 
systems, it is possible to steer any eigenstate of the free 
Hamiltonian to any neighborhood of any other eigenstate 
with arbitrary small norm of the control. 

5.2 Future works 

The estimates given for the norm only depend on the 
control potential B (and not on the eigenvalues of the 
free Hamiltonian A as long as the transition stay non- 
degenerate). It is possible that costs, with p < 1, are 
physically more relevant. A new approach would be needed 
to study this case since our methods do not provide any 
information about the norm of the control needed to 
steer an energy level to another in the case p < 1. 
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